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ABSTRACT 

We have obtained Magellan/IMACS and HST/ACS imaging data that resolve red giant branch stars 
in the stellar halo of the starburst galaxy NGC 253. The HST data cover a small area, and allow us 
to accurately interpret the ground-based data, which cover 30% of the halo to a distance of 30 kpc, 
allowing us to make detailed quantitative measurements of the global properties and structure of a 
stellar halo outside of the Local Group. The geometry of the halo is significantly fiattened in the same 
sense as the disk, with a projected axis ratio of b/a w 0.35 ± 0.1. The total stellar mass of the halo 
is estimated to be Mhaio ~ 2.5 ± 1.5 x 10^ Mq, or 6% of the total stellar mass of the galaxy, and 
has a projected radial dependence that follows a power law of index —2.8 ± 0.6, corresponding to a 
three-dimensional power law index of '~ —4. The total luminosity and profile shape that we measure 
for NGC 253 are somewhat larger and steeper than the equivalent values for the Milky Way and M31, 
but are well within the scatter of model predictions for the properties of stellar halos built up in a 
cosmological context. Structure within the halo is seen at a variety of scales: there is small kpc-scale 
density variation and a large shelf-like feature near the middle of the field. The techniques that have 
been developed will be essential for quantitatively comparing our upcoming larger sample of observed 
stellar halos to models of halo formation. 

Subject headings: galaxies: halos — galaxies: individual (NGC 253) — galaxies: structure — galaxies: 
formation — galaxies: spiral — galaxies: starburst 



1. INTRODUCTION 

In the favored hier archical model of galaxy formation 
(j White fc Reesl[T978| ). stellar halos are composed of the 
shredded remains of the smaller protogalac tic clumps 
that assembled the larger galaxy (Scarlc & ZinnI 119781 : 
iBuUock et al.|[200l|) . At radii approaching the virial ra- 
dius, the dynamical time within the halo approaches the 
age of the Universe. Stellar halos thus retain both spatial 
and kinematic information about this assembly history 
that is otherwise difficult to obtain. 

The stellar halo of our own galaxy has been extensively 
surveyed. Studies of high velocity stars in the solar 
neighborhood and bright stars at large distances have 
revealed a halo that follows a power law density dis- 
tribution, a total luminosity of approximately 10^ Lq, 
and abundant spatial, kinematic, abundance, and 
stellar population substructure (e.g. iHelmi et al.l IT999I: 



Maiewski et al.ll2003l:ryannv et al.ll2003HBelokurov et al l 



2007t ICarollo et alJ l2007t [Bell et al.l 120081 \20m 
Xue et al.l 120101: ICooper et al.i l2010aD . However, our 
position within the disk of the Milky Way makes it 
difficult to obtain a clear picture of the global structure 
of the halo. This is additionally complicated by the 
fact that different tracer populations must be used to 
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determine the structure at different distances. 

Some information has been obtained about the 
halos of external galaxies by searching for surface 
brigh tn ess features in d e ep ira ag ing (e.g. IBeck et al 
1982 IMalin fc Hadleyl \l99% IZheng et al.l 1199 
Martmez-Delgado et alT I2008L boiof T These obser- 
vations have identified prominent tidal streams that 
buttress the hierarchical model of galactic halo forma- 
tion. However, the fundamental limits in such studies 
are flat fielding errors and the brightness of the back- 
ground sky, which make global halo properties and very 
faint surface brightness features inaccessible (jde Jong! 
I2OO8I) . 

In order to overcome these difficulties, a new tech- 
nique has recently been employed: resolving individual 
bright stars, particularly red giant branch (RGB) stars 
which are present in all stellar populations, in the ha- 
los of external galaxies. Within the past few years, 
knowledge of the halo of M31 has blossomed with this 
technique, including the discovery of a metallicity gra- 
dient, rich spatial and kinematic substructure, and a 
global profile that is remarkably similar to that of the 



Milky Wa y ("Chapman et a l.l 120061: 



Ibata et al. 2007; McConnac hie et al 12009; Tanaka et al.l 
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2oToir Detections of halo stars in more distant galaxies 
are now beginning to be achieved, such as in individ- 
ual s mall HST/ ACS fields of NGC 5128 (IReikuba et al.l 
I2OO5I) and M 81 (IDurrefi et al.llMol) . and in NGC 891 us- 
ing ground-b ased Subaru/Suprime- Cam data over a large 
field of view (jMouhcine et al.ll2010D . 

One particularly interesting target is NGC 253. This 
nearly edge-on starburst galaxy, one of the largest galax- 
ies in the "Sculptor group" (which is not truly a bound 
group , but rather a fila.ment extended along the line of 
sight; IJerien et al.lll998l : IKarachentsev et aLll2003D . is of 
similar luminosity to the Milky Way and M31, allow- 
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ing direct comparisons. Deep optical images have long 
shown an optical halo extending beyond its disk, with 
evidence for a she l f- like feature to the so ut h of the disk 
([Beck et all ITosl IMalin fc HadlevI Hool . iFitzgibbonsI 
(|1990( ). in a photographic survey of the galaxy, also de- 
tected extended structure to the south of the galactic 
disk, particularly in the / band. He found a total light 
profile in the halo that declines as a power law with an 
index of ^ —2.5 and ~ —3.5 in the V and I bands re- 
spectively. 

Resolved near- infrared stellar photometry of t he south- 
east q uadrant of NGC 253 was per formed by iDavidgd 
see also iComeron et al.l 120011 who detected blue 
main sequence stars in the halo near the minor axis of 
the galaxy, most likely associated with the outflow from 
the central starburst). He detected an extended flattened 
distribution of asymptotic giant branch (AGB) stars ex- 
tending out to 13 kpc from the disk plane. Based on 
the similarity of the disk and halo -fC-band luminosity 
functions, which implies an extended star formation his- 
tory, he concluded that these stars were deposited into 
the halo via the disruption of the galactic disk in an in- 
teraction. The data were not of sufficient depth to detect 
RGB stars, which would be far more numerous and trace 
stellar populations of all ages. 

Theoretical models of galaxy formation in a ACDM 
universe have been used to predict the properties of stel- 
lar halos, usually with the assumption that the stars 
that are today found in the halo were accreted in the 
form of satellite galaxies that were subsequently tidally 
shredded (j BuUock et al.ll200 It [Bullock fc Johnston| [2005l: 
iPurceh et al., ,2007t iCooper et al.ll2010b[ ). These models 
predict that halos contain several percent of the total 
stellar mass of the galaxy, and have very clumpy struc- 
ture, particularly in the outer regions. They also predict 
significant galaxy-to-galaxy scatter in the total size and 
morphology of the halo structure (iBuUock fc JohnstonI 
[20051 IBeh et al.ll2008l: [Johnston et al.ll2008n . While these 
properties qualitatively match observed galactic halos, 
quantitative measurements of a larger sample of halos 
is required in order to adequately test these models and 
distinguish between them. 

Our goal is to test these models by obtaining resolved 
data of stellar halos in a larger, and therefore more rep- 
resentative, sample of galaxies. This will allow us to 
accurately determine stellar halo structural properties, 
and to investigate their correlation with the global pa- 
rameters of the associated galaxies. The larger sample 
will help determine any scatter between systems due to a 
stochastic formation history. Because of the large degree 
of predicted and observed substructure, obtaining an un- 
biased view of the halo requires a large field of view that 
is only possible using ground-based data. However, iden- 
tification and accurate photometry of stars in galaxies 
beyond the Local Group is far superior with HST. We 
are therefore using a combination of HST and ground- 
based wide-field imaging. 

In this paper, we present the first results of this pro- 
gramme, a study of NGC 253, along with a detailed de- 
scription of the analysis required to turn the data into 
quantitative global measures of halo structure. We adopt 
a distance modulus to NGC 253 of m - M = 27.71 {D = 
3.48 Mpc) throughout, as determined using GHOSTS 
(jRadburn-Smith ct al...2011i ). and consistent with but to- 



wards the upper end of distance estimates based on a va- 
riety of methods catalogued in the NASA Extragalactic 
DatabascQ (NED). Note that we use the term "halo" to 
refer to the stellar halo of the galaxy, not the dark matter 
halo, unless specifically stated otherwise. 

2. DATA 
2.1. Observations and Reduction 
2.1.1. Magellan/IMACS Data 

The ground-based data were obtained at the Magel- 
lan 1 (Baade) 6.5m telescope using the Inamori-Magellan 
Areal Camera fc Spectrograph (IMACS), a mosaic of 
8 2k X 4k CCDs, in f/2 configuration. This results 
in a O'.'2/pixel plate scale. The field was centered 14' 
south of the center of NGC 253, and therefore included 
a significant fraction of the galactic disk within the 24' 
unvignetted field of view. At the assumed distance of 
NGC 253, the field extends 30 kpc from the galaxy cen- 
ter. On the nights of 17 - 18 October 2009, we obtained 
1800 seconds in CTIO /-band and 2650 seconds in Bessel 
F-band. A further 1000 seconds were obtained in F-band 
on the night of 20 May 2010. The seeing fuU width hah 
max (FWHM) was approximately I'.'O - 1'.'2 during the 
exposures. Landolt photometric standards were taken on 
these and other nights during the same runs. 

Data were reduced using standard IRAF routinetlB, 
assisted by scripts developed for reducing Maryland- 
Magellan Tunable Filter data0. Note that although the 
IMACS documentation suggests using the column and 
row overscan regions instead of bias frames, we found 
that in the case of broadband imaging, saturated stars 
near the edges of the chips can bleed into the overscan 
regions, making them unusable. We therefore used bias 
frames that were obtained during the runs. Sky sub- 
traction was complicated by the presence of the galac- 
tic disk in the field, which made it impossible to fit an 
accurate low-order surface to the sky. We therefore de- 
termined the sky level using a two-dimensional median 
filter of width 101 pixels; this width was chosen to be 
large enough that photometry tests on individual stars 
found no bias, but small enough that the sky was flat 
to within 1% on large scales. This results in good image 
quality over most of the frame, but clearly invalidates the 
portion of the field that includes the disk of NGC 253, 
which contains smooth features on scales much larger 
than the filter width. However, as the goal of this work 
is to study the halo of the galaxy, this does not affect our 
results. An astrometric solution ove r the field wa s com- 
puted using the USNO A2.0 catalog (jMonet et al.i ri998'). 
A third order polynomial solution was required in order 
that the astrometric residuals did not vary systemati- 
cally across the field; the resulting solution achieved an 
rms scatter of 0'.'24, comparable to the pixel size. Two 
sample regions of the final stacked and sky-subtracted 
/-band image are shown in Figure [TJ one located ~ 9.'5 
(10 kpc) from the center of the galaxy and one located 
at ^ 22' (22 kpc), both approximately along the minor 

http:/ /nedwww. ipac.caltech.edu/ 

* IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 

^ http:/ /www. astro. umd.edu/~veilleux/mmtf/datared. html 
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Figure 1. Left: Sky-subtracted stacked IMACS /-band image of 
a 82" X 110" region located ~ 9.'5 (10 kpc) from the center of 
NGC 253, along the minor axis. RGB stars within its stellar halo 
are clearly visible as faint point sources. Right: As in the left panel, 
for a region located ~ 22' (22 kpc) from the center of NGC 253. 
Very few if any RGB candidates are visible. 

axis. The faint point sources that are abundant in the 
inner field but almost completely absent in the outer field 
are RGB stars in the halo of NGC 253. We ignore the 
very outermost part of the mosaic, which is compromised 
by the shadow of the guide star camera, whose location 
varied between exposures. 

2.1.2. GHOSTS 

We use data f r om the HST/GHOSTS survey 
(|de Jong et all 120071: IRadburn-Smith et al.l 120 111) over- 
lapping with our IMACS pointing for the purposes of car- 
rying out our artificial star tests, and to test our inferred 
stellar densities (in the high stellar density regime). We 
use HST/ACS F606W and F814W data for N GC 253's 
Fields 9 and 10 (|Radburn- Smith et al.ll2011[) . with ex- 
posure times of 2 X 680s in both F606W and F814W 
bands. Data were reduced and calibrated using the 
standard HST data reduction, and photometry was per- 
formed using the ACS module of DOLPHOTO. Photo- 
metric parameters were converted to V an d /-band mag- 
nitude s using the conversions presented in lSirianni et al.l 
(|2005f ). Typical uncertainties (random and systematic) 
are < 0.05 mag at the / < 26 limits of interest in this 
paper. 

2.2. Photometry 

Sources were found on the stacked IMACS /-band im- 
age using the IRAF DAOFIND task, with a threshold 
of 3 times the typical standard deviation of the back- 
ground. Aperture photometry was performed using the 
IRAF APPHOT package on the stacked V- and /-band 
images at the locations found by DAOFIND, using an 
aperture of radius 8 and 6 pixels respectively, matched 
to the width of the point spread function (PSF) to max- 
imize the signal-to-noise. Aperture corrections were cal- 
culated in each band using 5 bright unsaturated rela- 
tively isolated stars. Magnitudes were de-reddened using 
the Galactic foreground extinction values catalogued in 
NED for NGC 253. The final catalogue consists of all 
stars that have a positive measured fiux in both bands. 

Photometric calibration was performed using aperture 
photometry of standard stars in the Landolt fields SA98, 
SAlOl, and SA114 taken on the nights of 17 - 19 October 

http:/ /purccU. as. arizona.edu/dolphot/ 




20OO 4000 6000 8000 

X 

Figure 2. The locations of all stellar sources are plotted in pixel 
coordinates on the bottom-left axes and angular coordinates on the 
sky on the upper-right axes. Points within the circle are considered 
far enough from the edge of the field that they are usable (usable 
points are colored red in the online version of this figure). Out- 
lines of the footprint of the IMACS chips for one particular dither 
position are shown by the rectangles. 

2009, using an aperture of 40 pixels, sufficient to enclose 
the entire fiux of each star while not overlapping with 
neighboring stars. We fit a zero point, as well as airmass, 
color and radial term^ to the residuals. The airmass 
term in the /-band and the color terms in both bands 
were found to be negligible, and so were assumed to be 
zero. The resulting photometry had an rms dispersion of 
0.051 mag in /-band and 0.054 mag in V-band, constant 
over the 10 < V < 15.3 magnitude range of the standard 
stars. As will be seen later, such a level of photometric 
calibration accuracy is sufficient for our purposes, as it 
is smaller than other sources of systematic uncertainty. 

A plot of the locations of all sources is shown in Fig- 
ure [21 The points within the circle were used; this circu- 
lar boundary is also used in the Monte Carlo simulations 
discussed in Section 13.21 to determine the effective sur- 
veyed area. Although the dithering pattern fully covered 
the chip gaps, the effective exposure time in the areas 
covered by chip gaps is lower and therefore the noise 
level is higher. This is reflected in the larger number 
of counts above the threshold in those regions, since the 
fixed threshold is no longer equal to 3cr. However, the 
excess counts are predominantly at the faint end, fainter 
than the brightness of the RGB candidates. The impact 
of the chip gaps is therefore small, although it is large 
enough that artificial structure along the horizontal chip 
gap is faintly visible in maps of RGB stars. 

3. COLOR-MAGNITUDE DIAGRAMS AND ARTIFICIAL 
STAR TESTS 

In this section, we present the color-magnitude dia- 
grams (CMDs) of the ground-based and HST/ACS data, 
identify RGB stars in the data, and use artificial star 
tests to understand the relationship between the RGB 
star counts derived using the ground-based data and the 
intrinsic stellar population. 

3.1. Color- Magnitude Diagrams 

In Figure[5^, we present the CMD of all recovered point 
sources within the main IMACS field. The bright blue 

The radial term accounts for the increase in solid angle sub- 
tended by pixels at larger radii. 
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Figure 3. Color-magnitude diagrams of the Magellan/IMACS field (a) and GHOSTS HST/ACS field 10 (b) around NGC 253. The 
absolute magnitude scale on the right-hand axis assumes a distance modulus of m — M = 27.71. 



vertical feature is primarily composed of pieces of the 
main body of NGC 253 that have been selected as stars 
by DAOFIND, along with the shredded diffraction pat- 
tern of foreground Milky Way main sequence stars, and 
background disk galaxies. At the bottom, a wide feature 
that appears to be the RGB is seen. However comparison 
with the GHOSTS Field 10 data (Figure[33) reveals that 
this putative RGB extends to much brighter magnitudes 
than in the HST/ACS data. Comparison of the IMACS 
and HST catalogs over their area of overlap reveals that 
the majority of "stars" detected in the IMACS catalog 
are blends of two or more stars. Therefore, while the 
detections in IMACS are indeed detections of RGB stars 
in the halo of NGC 253, interpretation of the number 
counts requires artificial star tests in order to determine 
the relationship between the numbers and properties of 
stars in the IMACS catalog and the true numbers and 
properties of the stellar population. 

3.2. Artificial Star Tests 

In order to test the photometric accuracy and recov- 
ery/contamination rates of TRGB stars in NGC 253, we 
carried out a suite of artificial star tests. Artificial stars 
with mi < 26 were placed into a Ik x Ik region in the 
outer parts of the V and /-band NGC 253 images with 
empirical PSFs generated from well-exposed stars. V- 
band and /-band apparent magnitudes were Monte-Carlo 
sampled from the CMD of the HST/GHOSTS NGC 253 
Field 10 pointing. The stars were distributed randomly 
across the region at 8 different densities (labelled "30k" 
through "4000k" depending on the number of artificial 
stars that were added), spanning the observed range of 
RGB star densities seen in the IMACS image and in the 
GHOSTS data. The images were then output and pro- 
cessed by the photometric pipeline in the same way as 
the data. 

Photometry was performed on the artificial star frames 
exactly as for the original frames. CMDs generated from 
the input magnitudes and from the output photometry 
on the artificial star frames for three sample artificial star 
tests (the 30k, 700k and 4000k tests) are shown in Fig- 
ure m As the stellar density increases, the RGB stars 
become blended and are detected at magnitudes signif- 
icantly brighter than those of individual stars. This re- 
sults in a CMD exactly like the one observed if the ma- 
jority of our sources are found in regions of high density. 



In order to determine the number density of stars at 
the tip of the red giant branch (TRGB), we have com- 
pared the number of input RGB stars above a threshold 
of / < 24.25 (corresponding to an absolute magnitude of 
Mi < -3.46 at the distance of NGC 253) to the num- 
ber density of output detections within a box extending 
from / = 23 to 24, and from - / = 1.0 to 2.5. This box 
lies entirely above the input threshold in order to better 
capture the bulk of the halo stars at the typical densi- 
ties of our data (see Figure S]) . We have tested several 
different criteria for both the TRGB threshold and the 
magnitude boundaries of the selection box, ranging by 
±0.5 mag. We find that this introduces ^ 20% changes 
in the derived number densities once they are normal- 
ized by the expected number of stars above each TRGB 
threshold for an old metal-poor population, as discussed 
in Section [5l 

The relationship between the input number density 
above the TRGB threshold and the output number den- 
sity within the selection box is shown as the right-hand 
(high-density) dashed line in the upper panel of Figure [SJ 
while their ratio (i.e. the recovered fraction) is shown 
in the bottom panel. A problem is immediately appar- 
ent: because there are 183 detections within the RGB 
selection box in the "blank field" to which the artificial 
stars were added (a combination of background galax- 
ies and real RGB stars in the stellar halo), more stars 
are detected than were added at low densities. This is 
a property of the baseline image, not of our ability to 
detect stars, and so is artificial. One could simply sub- 
tract 183 from the number of output stars (shown as 
the left-handed dashed line in Figure [S|), but this is also 
incorrect: at high densities, we approach the confusion 
limited regime, and the original detections have mostly 
been covered up by artificial stars and are therefore no 
longer detected. A better approach is to compare the 
locations of the detections in the artificial star frames to 
the locations of the original detections in the blank field. 
We performed a globally-optimized one-to-one match of 
the closest original detection to each detection in the ar- 
tificial star frame, and omitted any whose best match 
was within 0'.'5. This results in the black solid line in 
Figure [S] The fraction of blank field sourcees that are 
omitted, /bg.omit, is a useful measure of the degree to 
which individual sources are likely to be covered up. 

At the confusion limit, the number of output detec- 
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tions saturates, and even begins to decrease as the ma- 
jority of the blends become brighter than the selection 
box. We therefore set a threshold of lognout = 5.60, 
at which we can only set a lower limit on the true den- 
sity (logriin > 6.39), shown as the right-hand vertical 
gray line in Figure El Fortunately, very little of the data 
reach this limit. At the low density end, we note that 
the recovered fraction approaches unity, and so we make 
the ansatz that the recovered fraction has the form of a 
decaying exponential at densities below that of our 30k 
artificial star test; this results in the gray line extending 
to the left, which we use to extrapolate to lower densi- 
ties. Again, most of our data lie above the point where 
we begin using this extrapolation, at lognout < 4.39 
(logriin < 4.46). Our final correction from the artificial 
star test is therefore the combination of the solid black 
and gray lines in the top panel of Figure [5l Interpolation 
is performed linearly between logriout and lognin. The 
error bars in the bottom panel indicate the 20% uncer- 
tainty due to the choice of TRGB threshold magnitude 
and CMD selection box. 

Note that in order to convert number counts to densi- 
ties, we required the effective area of the observed region. 
This is also required in order to calculate the densities 
in annuli, sectors, pixels in a density map, and the foot- 
prints of the GHOSTS pointings. Whenever the effective 
area was required, Monte Carlo simulations were per- 
formed. A^MCtot = 10^ points were randomly sampled 
uniformly over a patch of the unit sphere of solid an- 
gle fiMCtot = 0.2088n° encompassing the entire IMACS 
field. Within any given region, we counted the number 
of Monte Carlo points lying within that also lie within 
the field where stars were chosen (i.e. within the circle 
on Figure [2]), A'mCi and also the number of bad pixels, 
A'bad- The bad pixels included not only pixels formally 
counted as bad during the data reduction, but also re- 
gions that were completely devoid of detections due to 
bleed trails or the presence of a bright foreground star. 
We assume that the effective solid angle of each bad pixel 
is rJbad = 0.04n", as given by the formal IMACS f/2 plate 
scale; the true effective area varies over the field, but this 
effect is small relative to the two order of magnitude vari- 
ation in the number density of RGB stars across the field, 
particularly given that the fraction of bad pixels is small. 
The total solid angle of the region is then given by 

i^rcgion = A^MCTf— — — A^bad^^bad- (1) 
jVMCtot 

3.3. Background Galaxy Subtraction 

Some fraction of the detections must be due to back- 
ground galaxies. To assess the impact of this back- 
ground, we have used a control field that was created 
using V and /-band imaging data taken using the Wide 
Field Imager on the ESO/MPG La Silla 2.2-m tele- 
scope of the Extend ed Chandra Deep Field South (see 
iGawiser et al.l I20060 . These calibrated images cover a 
slightly larger field of view than the IMACS data, and 
are somewhat deeper and have a smaller point spread 
function (0'.'9). The data were resampled onto 0'.'2 pixels 
(to match IMACS), convolved with the correct Gaussian 
to bring the FWHM to 1'.'04 and 1'.'2 for 7-band and V- 
band respectively, rescaled to the same zero point as the 
IMACS data, and sky noise was added to bring the RMS 



of the image into agreement with the IMACS data for 
NGC 253. We then ran the images through our pipeline 
for star detection and photometry, and used the output 
photometric catalog as our "control field" . The density 
of sources within the RGB selection box in the control 
field is logncontioi =4.27n°~^. We scale this density by 
/bg,omit before subtracting it; at densities below that of 
the 30k artificial star frame we subtract the full control 
density. The background counts are much smaller than 
the densities we derive over much of the halo, but are 
comparable to the densities in the outermost regions of 
the field. 

4. HALO STRUCTURE 
4.1. Profiles 
4.1.1. Radtal Profile 

The IMACS field was divided into radial bins of width 
4 kpc, excluding the region within 5 kpc of the disk plane 
in order to avoid thin and thick disk stars. These annuli 
were further subdivided into angular sectors of azimuthal 
extents /S.9 — 20°, in order to determine the structure of 
the halo at a given radius (different azimuthal bin sizes 
ranging from 16° to 40° were also tested, with no signif- 
icant differences in our conclusions). The top- right inset 
in Figure [HI shows the sectors; azimuth is defined with 
= 0° along the major axis to the north-east (upper-left 
in the Figure) and increasing counter-clockwise. 

The radial projected density profile is shown in the 
main panel of Figure |6l The error bars are the sum, in 
quadrature, of the Poisson error in the number of de- 
tected stars in each bin and the 20% uncertainty due 
to our choice of TRGB threshold and CMD selection 
box (see Section [32]) • The profile is well fit by a broken 
power-law, of inner slope n = —1.0 ± 0.4 and outer slope 
n — —4.0 ± 0.8 (where the quoted error is the formal un- 
certainty in the fit, which is an underestimate; see below), 
with a break at 14.9 kpc. As demonstrated later, the 
break between these power-laws is not due to an intrin- 
sic change in the distribution of halo stars, but is rather 
an artifact of using circular annuli to describe a distri- 
bution that is intrinsically flattened. The IMACS counts 
within the GHOSTS fields are shown as the triangles, 
while the HST/ACS counts are shown as the squares. 
The difference between the triangles and squares gives 
an indication of the systematic uncertainty, ^ 0.2 dex, 
which dominates over the statistical errors shown in the 
profiles. This is discussed in more detail in Section 14.31 
Given the size of this systematic uncertainty, a more re- 
alistic error for the slope of the radial power law can be 
obtained by shifting the inner points up by 0.2 dex; doing 
this changes the slope of the inner power law to n = —1.6 
and the the slope of the outer power law to n = —5.3. 
Although it appears that the IMACS counts within the 
outer GHOSTS field (Field 10) are much higher than the 
counts at adjacent radii, this is because the radial profile 
averages over all azimuths, while the GHOSTS field is at 
one specific azimuth. Indeed, as shown by the azimuthal 
profiles in the bottom-left inset of Figure [6] (also shown 
as the gray lines with error bars on the main plot; these 
gray points are offset horizontally with azimuth for clar- 
ity), there is a strong trend in density with azimuth for 
each radial bin. The IMACS counts in the GHOSTS field 
are fully consistent with the counts in the same radial 
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Figure 4. Artificial star CMDs. The left column shows input magnitudes of the artificial stars that were placed in the artificial star 
frames, while the right column shows the CMDs determined by performing photometry on the output artificial star frames. In the output 
CMDs, red points are those located within O'.'S of a star in the original image. The rows correspond to injecting 30k (top), 700k (middle), 
and 4000k (bottom) artificial stars. The input and output density of stars is given in each panel, in number per square degree. The lines 
in the left column and the boxes in the right column denote the TRGB magnitude threshold and the boundaries of the RGB selection box 
respectively. 

bin when this trend is taken into account. OveraU, the 
almost uniform positive gradient in the azimuthal pro- 
files indicate that the counts towards the galactic plane 
(at larger azimuth, which is defined to increase counter- 
clockwise) are significantly higher than along the galactic 
minor axis at the same radius. This is strong evidence 
that the halo profile is not spherically-symmetric, but is 
rather flattened in the direction of the disk. Another 
possibility is that we are not seeing a true stellar halo, 
but rather a very extended thick disk, which would nec- 
essarily have a very flattened distribution (see discussion 
in Section H. 1 .31 below) . 
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Figure 5. Results of the artificial star tests. The x-axis gives 
the density of stars recovered within the RGB selection box in the 
artificial star frame, while the y-axis gives the density of artificial 
stars injected into the artificial star frames that lie above the TRGB 
threshold (top), and the ratio of output to input stars (bottom). 
The dotted line in the upper panel denotes riin = riout. The right 
dashed line shows the case where no correction is made for the 
presence of stars in the baseline image. Subtracting the density of 
stars found in the original image gives the left dashed line, while 
subtracting only the stars in the original image that are located 
within 0'.'5 of a star in the artificial star frame gives the adopted 
solid line. The gray lines show the adopted extrapolations. The 
error bars in the bottom panel indicate the 20% uncertainty due to 
the choice of a particular TRGB magnitude threshold and CMD 
selection box. 



4.1.2. Elliptical Radial Profile 



The azimuthal trend of the halo density in Section l4.1.1l 
strongly indicates that the stellar halo of NGC 253 is 
flattened. Therefore, in order to determine an accurate 
radial profile, we must adopt "annuli" that mimic the 
intrinsic distribution of the stars. In this section, we 
adopt elliptical annuli aligned with the stellar disk. For 
the elliptical axis ratio, we have tested values r anging 
from the observed disk axis ratio of 6/a = 0.25 (jPencd 
[19801 iKoribalski et al.|[200l . up to b/a = 0.5. In the 
azimuthal direction, we adopt an angular coordinate (j) 
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Figure 6. Main panel: The radial profile of the proj ected TRGB 
stellar density in the IMACS field, as defined in Section [3.2l The di- 
amonds indicate the average value within each radial annulus. The 
dark solid line is a double power-law fit to the radial profile, which 
has a break at 14.9 kpc. The continuation of the outer power-law 
to smaller radii is shown as the dotted line. The TRGB number 
densities found in the HST/ACS GHOSTS data are denoted by the 
squares, while the densities found in the same fields using IMACS 
data are shown as the triangles. The azimuthal profile within a 
given radial annulus is shown as the gray line at each radius. The 
plotted error bars for both profiles account for the Poisson error 
and the uncertainty due to our choice of TRGB threshold and CMD 
selection box; however, the true uncertainty is dominated by sys- 
tematics, the magnitude of which can be estimated by comparing 
the squares and triangles. Top-right inset: The points indicate the 
locations of RGB sources in the IMACS field, projected onto the 
sky at the distance of NGC 253. The coordinate system is labeled 
in kpc from the center of the galaxy, with north up and east to 
the left. The field is divided into radial annuli of width 4 kpc, and 
further subdivided into sectors with azimuthal extent A8 = 20°, 
as indicated by the coloring. The squares indicate the locations of 
the GHOSTS fields. Bottom-left inset: Azimuthal profiles of the 
projected TRGB stellar density. The x axis gives the azimuthal 
angle, 8, in degrees. Each line corresponds to a different radius, 
starting at the top and going down. The points are colored the 
same as the corresponding sectors in the top-right inset. The tri- 
angles and squares are as in the main panel. Azimuth increases 
counter-clockwise in the top-right inset, so the right-most angular 
sector in this panel is the one closest to the disk plane. 

defined by 

tan0=^-, (2) 

where rmaj and are the projections of the radial vec- 
tor along the major and minor axes respectively. If the 
halo distribution were perfectly disk-like, angular sec- 
tors spaced uniformly in (j) would correspond to equal 
volumes. Our results are qualitatively unchanged if we 
instead adopt sectors that are spaced uniformly in the 
true projected azimuth 9. The elliptical annuli had mi- 
nor axes spaced by 4 kpc, which were further separated 
into sectors with azimuthal extent A0 = 8°, as shown 
in Figure [7] (A0 values from 4° to 10° showed similar 
behaviour). 

The radial halo profile in elliptical annuli with axis ra- 
tios b/a = 0.35 is shown in the main panel of Figure [71 
A single power-law with a slope of —2.8 ± 0.3 provides 
an acceptable fit to the data. This further demonstrates 
that the break in the power law seen in Figure [6] is en- 
tirely an artifact of the assumed geometry: the circular 
annuli coupled with the 5 kpc inner cutoff probe different 
ranges of elliptical annuli as a function of radius. Shift- 
ing the inner points up by 0.2 dex, as in Section [4.1.11 
steepens the slope lo n = —3.4. 

The azimuthal profiles are shown in the bottom-left in- 
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Figure 7. As in Figure (6] but using elliptical annuli of axis ratio 
b/a = 0.35, with minor axis widths of 4 kpc and elliptical sectors 
with azimuthal extents of = 8° . A single power law is fit to the 
radial profile. The x axis of the bottom-left inset gives the angular 
coordinate, </>, in degrees. 

set of Figure [7| Unlike with the circular annuli, the den- 
sities do not show a systematic gradient, again demon- 
strating the flattening of the halo. However, the az- 
imuthal profiles are far from flat, and exhibit a strong 
peak at intermediate azimuths, (f) ^ 100°. This could 
either indicate that the halo isophotes are smooth but 
"boxy" , or could be an indication of discrete structure 
within the halo. The fact that the peak in azimuth is 
broad, rather than being confined to a single bin, indi- 
cates that it is not an artifact of the chip gap (which 
runs through many of the bins where the density peaks, 
particularly at large radius) or a manifestation of a thin 
tidal stream; it must be due to an extended or smooth 
fe ature. The shelf seen i n the deep optical photometry 
of [Malin fc HadlevI ([IQQTT ) extends in this direction, and 
it is likely that our azimuthal peak corresponds to this 
structure. 

The choice of axis ratio b/a — 0.35 was chosen by eye to 
minimize the variation of the azimuthal profiles. Values 
ranging from the disk el lipticity of b/a = 0.25 to the 
value of b/a — 0.4 that [Davided (|2010D estimates for 
the stellar halo all provide reasonable, though somewhat 
poorer, descriptions. 

4.1.3. Vertical Profile 

In the previous section, we found that the geometry of 
the halo must be quite flattened, being well fit on ellipses 
that are almost straight lines parallel to the disk over 
much of the field. It is therefore natural to ask whether 
a pure disk description, with the density depending only 
on distance from the galactic plane, provides a better de- 
scription. We have constructed a vertical profile, calcu- 
lated the density within vertical slices of Az — 4 kpc and 
further subdivided into rectangles of length 4 kpc in the 
direction parallel to the major axis, as shown in the top- 
right inset of Figure [51 (major axis subdivisions ranging 
from 2 to 8 kpc gave similar results). The innermost ver- 
tical sheet is somewhat compromised by incompleteness 
within the disk near the minor axis, and so the density 
of this bin was not used when fitting the density profile. 
The shapes of the ellipses and of the vertically-stratified 
sheets are not too dissimilar over much of the field, and 
so it is not surprising that both the quality and slope 
(—2.8 ± 0.2, steepening to —3.3 if the inner points are 
systematically low) of the best fit power law are similar 
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to the elliptical case. The density variation parallel to 
the major axis within each layer is also similar to the 
angular variation within the elliptical annuli. 

Given that a description where the density varies pri- 
marily in the direction perpendicular to the disk is ac- 
ceptable, it is natural to wonder whether the extended 
stellar structure we have detected is truly the stellar halo, 
or is rather a thick disk. This is pertinent given that 
M31, the galaxy whose outer regions are the best studied, 
shows evide nce for both an e xtended rotating disk- like 
component (jibata et al.l 12005 1 ) and a thick d isk of esti- 
mated scale height 2.8 kpc (jcijllins et ani2011f ). In order 
to test this hypothesis, we have fitted a vertical expo- 
nential to the density structure we have detected (again 
omitting the innermost point), shown as the dashed line 
in Figure [51 An exponential is found to fit the data as 
well as the power law, but with an extremely large scale 
height of hz = 5.3 ± 0.5 kpc (4.5 kpc if the inner points 
are moved systematically higher). While we cannot rule 
out this possibility with the present data, we believe that 
a power law halo is a more likely explanation for the fol- 
lowing reasons: 

1. Such a large scale height is significantly larger than 
that of other stellar disks of comparable galaxies. 
While it is only 1.9 times larger tha n the 2.8 kpc 
thick disk scale height estimated bv iCollins et al.l 
()2011D for M31, these authors did not directly mea- 
sure the scale height but rather assumed that the 
scale height varies with scale length in the same 
manner as the sa mple of low surface brightnes s 
galaxies studied bv lYoachim fc DalcantonI ()2006D . 
The large thick disk scale height inferred for M31 
is thus a direct consequence of its long disk scale 
length. NGC 253, on the other hand, has a 
much smaller scale length; estimates based on sur- 
face photometry at wavelengths from the opti- 
cal through far infrared place the thin disk scale 
length between 1.5 and 2.6 kpc, once sca led to 
our adopted distance of D = 3.48 Mpc (IPencd 
[19801: IForbes fc Depovl[l99llRadovich et al.ll200lli 7 
Adopting a thin-to-thick disk scale length ratio 
of 1.3 and the thic k disk scale heig ht-to-scale 
length ratio of 0.35 (iCollins et al.ll201lL based on 
lYoachim fc DalcantonI I2006D . the expected thick 
disk scale height is between 0.7 and 1.2 kpc, many 
times smaller than what we derive. The derived 
scale height would require a very large vertical ve- 
locity dispersion to maintain; for a stellar mass sur- 
face density in the disk of 75 Mp^ pc~^ (as for the 
Milky Way at the solar circle, iBinnev fc Tremaind 
[2001 this is a reasonable assumption given the sim- 
ilar total luminosity and disk scale length of the two 
galaxies), cr^ — 104 km s~^, making it essentially a 
pressure-supported structure. 

2. Both IMACS (GHOSTS fields) data points lie 
above the extrapolation of the exponential fit. 

3. The apparent quality of the exponential fit comes 
mainly from the innermost point, but this point is 
likely an underestimate due to incompleteness. 

This issue can be resolved with more data that are 
deeper and cover a larger area; new HST GHOSTS data 
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Figure 8. As in Figures [6] and [T] for bins of width 4 kpc perpen- 
dicular to the disk and 4 kpc along the disk major axis. The x-axis 
in the lower-left inset indicates the distance along the major axis 
from the center of the galaxy. The solid line is a power law fit to 
all but the innermost point, which is somewhat compromised by 
incompleteness within the disk. The dashed line is an exponential 
fit to the same data points, with scale height = 5.3 kpc. 

on NGG 253 and other galaxies are being analyzed to ad- 
dress the outer shape and profile of extended halo com- 
ponents. 

4.2. Density Map 

A map of the derived density of TRGB stars in the halo 
of NGC 253 is presented in Figure [5] This map was pro- 
duced by binning the stars into pixels 1 kpc on a side, 
and performing the background subtraction and artifi- 
cial star correction. The image was then stretched using 
histogram equalization to enhance visibility, smoothed 
using the minimum curvature surface algorithm (FranQ 
11982'). and resampled to a resolution of 0.1 kpc. The 
main body of the galaxy, where we cannot reliably de- 
tect stars, is apparent as the dark ellipse at the top (an 
optical image of the galaxy is overlayed, for reference), 
while the GHOSTS pointings are denoted by the black 
squares. Magenta indicates regions where the derived 
density reaches the saturation level in the IMACS data 
and we are therefore only able to place a lower limit on 
the true stellar density at these points. 

A few features are notable in this map. First and fore- 
most, the stellar density falls off strongly with radius 
from the galaxy, demonstrating that the stars belong to 
NGC 253. Secondly, the isopleths of constant stellar den- 
sity are ellipses whose axis ratio is similar to or somewhat 
larger than that of the observed galaxy, as anticipated 
from Section 14.1.21 Thirdly, there are are some large- 
scale substructures, chiefly a shelf-like feature near the 
center of the field (shown schematically by the dashed 
curve; see also Malin fc Hadlcy 1997 figure 4). Finally, 
there is a small scale pixel-to-pixel variation in the den- 
sity. None of these features appears to be associated with 
any overdensity of resolved background galaxies. Note 
that the vertical feature at a; = extending from y = —20 
to —25 kpc lies right along the chip gap (shown with the 
dotted line), and is therefore most likely an artifact. 

In Figure IIOI we show the residuals between the den- 
sity map and the best fit power laws, expressed as the 
logarithm of the ratio. The left-hand panel shows the 
result using the radial power law. The systematic trend 
for the residuals to be negative on one side and positive 
on the other side demonstrates again that a simple ra- 
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Figure 9. Map of the density of TRGB stars in the halo of NGC 253. The map is oriented with north up and east to the left. The x and 
y axes are labeled in kpc from the center of the galaxy. Magenta indicates where the derived density is a lower limit on the true density. 
The GHOSTS pointings are denoted with the black squares, the location of the shelf is shown by dashed lines, and the vertical dotted line 
denotes the horizontal chip gap, where the counts are artificially inflated. The dark ellipse at the top is the main body of the galaxy, where 
individual sources can no longer be distinguished; an optical image of the galaxy, to scale, has been overlayed for reference (image courtesy 
of T.A. Rector/University of Alaska Anchorage, T. Abbott and NOAO/AURA/NSF). The map was generated using 1 kpc X 1 kpc pixels 
and then smoothed using the minimum curvature surface. Note that the mapping between the color scale and density value has been 
chosen to enhance visibility, so attention should be paid to the labels on the color bar, which give the logarithm of the TRGB number 
density per square degree. 
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Figure 10. Ratio between the density map shown in Figure|9]and the best fit radial (left), elliptical (middle) and vertical (right) power 
law from Figures [G] [3 and [8] respectively. The color scale label is the logarithm of the ratio. Axes are labeled in kpc from the center of the 
galaxy. 



dial power law is a poor description of the halo, i.e. it 
is not spherically symmetric. In the middle panel, the 
elliptical power law fit is used. In this case, the residuals 
show no systematic trend with azimuthal angle, demon- 
strating that the halo is consistent with being flattened 
in the same sense as the disk. The vertical feature at 
a:; = 0, which is an artifact of the chip gap, is very ap- 
parent, along with smaller-scale variation. In the right- 
hand panel, the vertical power law fit is used; given the 



similarity between the elliptical and vertical geometry, it 
is not surprising that this panel shows similar features 
to the middle panel. However, there is a small system- 
atic negative trend at the on the right side of the image, 
which again suggests that a power law halo is a better 
description than an exponential thick disk. 

4.3. GHOSTS Fields 

The overlap between the GHOSTS fields and the 
IMACS data provides a check on our methodology. The 
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density of TRGB stars brighter than / < 24.25 in the 
HST/ACS GHOSTS data is logTVTRGB = 5.96 and 
logA^TRGB = 6.13 for Fields 9 and 10 respectively, 
while the derived densities from the IMACS data over 
the identical areas of sky are logTVxRGB = 5.73 and 
logA^TRGB = 5.92, factors of 1.7 and 1.6 lower, respec- 
tively. While these factors are not unity, they are re- 
assuringly close given the dramatic difference in image 
quality and depth between the HST/ACS and IMACS 
data, and validates the use of moderate ground-based 
data to study stellar halos of galaxies beyond the Local 
Group. 

The factor of^^ 1.7 (~ 0.2 dex) provides an indication 
of the magnitude of uncertainty that is introduced by us- 
ing ground-based data. Because the factor is similar in 
both fields, it is not a random error but must reflect a 
systematic uncertainty, and dominates over the random 
Poisson error (which is denoted by the error bars of Fig- 
ures[n]-|S]). With only two overlapping HST fields, which 
are both located at similar densities, we cannot deter- 
mine whether this systematic reflects an overall normal- 
ization error, or is specific to the highest density regions. 
Overlapping HST data at lower stellar densities will be 
required to disentangle these possible effects. Although 
there is a small magenta region in Figure [S] within the 
GHOSTS Field 10 pointing, indicating that the IMACS 
data give only a lower limit at that location, the lack of 
magenta within the Field 9 pointing indicates that this 
is not the cause of the offset. 

Note that although the photometric accuracy of the 
IMACS data is insuffic ient to determine met a llicity via 
the color of the RGB, iRadburn-Smith et al.l (|201lb re- 
port an [Fe/H] estimate based on the TRGB colors of 
-1.1 ± 0.3 for the HST/ACS GHOSTS data in both of 
these fields. 

4.4. Shelf 

As noted in Section an overdensity of stars is seen 
on the south side of the galaxy. A rough outline of the 
region is shown with the dashed curve in Figure 1^1 This 
overdensi ty coincides with the op tical shelf seen in deep 
images of IMalin &: HadlevI ()1997h . Our detection of this 
feature demonstrates that it is due to stars, rather tha n 
synchrotron emission (as suggested bv lBeck et al.lll982t ). 

Knowledge of the stellar content in the shelf would 
greatly aid in understanding its origin. Unfortunately, 
the high density in this region results in significant blend- 
ing in the IMACS data, and therefore errors in the stellar 
colors that are too large to allow the detection of stel- 
lar population differences. However, the locations of the 
GHOSTS pointings are fortuitous: Field 10 lies on the 
shelf, while Field 9 lies just off the shelf. These fields are 
at almost identical elliptical radii, so we would expect 
the stellar populations in these fields to be very similar 
if the halo were smooth. We can therefore compare the 
colors of stars detected in the HST/ACS data in these 
fields to constrain the stellar population of the shelf. 

We have plotted the fraction of TRGB stars (/ < 
24.25) as a function of their V ~ I color for the two 
GHOSTS fields in Figure HH The distribution is very 
similar, indicating that the stellar populations in the 
shelf are not dramatically different from the rest of the 
halo at this distance. The differences that are seen are in 
the sense that the shelf contains more red (1.8 < V — I < 
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Figure 11. Histograms of thie V — I color distribution of stars 
at the TRGB in GHOSTS Fields 9 ("OFF Shelf") and 10 ("ON 
Shelf"). 

2.6) TRGB stars, and fewer blue (1.2 < V - I < 1.8) 
TRGB stars. Using the tr ansformation between m etal- 
hcity and TRGB color in iBellazzini et al.l ()200H ). this 
corresponds to more stars at —1.3 < [Fe/H] < —0.7 

and fewer at [Fe/H] < -1 .3. Indeed, using the same 

data. ' Radburn-Smith et al.l (|2011| ) find that the Field 10 
TRGB is 0.06 mag redder and 0.05 dex more metal-rich 
than Field 9. 

A higher metallicity in the shelf is consistent with hier- 
archical models of stellar halo formation, which predict 
that the highest surface brightness substructures, having 
been formed from the largest satellites, are more metal- 
rich t han the well-mixed smooth component (jFont et al.l 
12001) . 

5. HALO LUMINOSITY AND MASS 

5.1. Synthetic CMD Calibration 

In order to translate the number of TRGB stars to 
a total luminosity of the unde rlyi ng stellar population, 
we used the online lAC- STAlQ synthetic CMD code 
(|Aparicio fc Gallartll2004D to generate an old metal-poor 
stellar population, and then compared the number of 
TRGB stars to the total luminosity. We tested sev- 
eral sets of parameters for the stellar population: a 
single burst lasting for 1 Gyr that ended 12 Gyr or 
9 Gyr ago, and a metalhcity oi Z = 0.0007 w 1/30 Zq, 
Z = 0.002 » 1/10 ZfTi, or Z = .007 « 1/3 Zq, in aU 
cases using a 'Kroupa et al.l ([IMS) initial mass function, 
the Girardi et al. (200(J) stellar evol ution library, and 
bolometric corrections based on the iCastelli fc Kuruc3 
(2004) models. The synthetic CMDs contained between 
10000 and 20000 stars above M/ < -3.46 (equivalent to 
our TRGB cut at an observed magnitude of / < 24.25). 
The resulting stellar populations had a total integrated 
luminosity per TRGB star of between 1.2 x 10^ and 
1.5 X lO'* Lq/Ntrgb, increasingly monotonically with 
metallicity but having very little dependence on the age 
of the population (L/A^trgb is ~ 3% higher for the older 
population). We adopt a scaling of 1.3 x 10** Lq/Ntrgb, 
appropriate for an old population with metallicity Z = 
0.1 Zq, close to the inferred metal licity in the overlap- 
ping GHOSTS fields padburn-Smith et al.ll201lD . and 
caution that large differences between the assumed and 
intrinsic metallicities could introduce ~ 20% systematic 

■"-^ http;/ /iac-star. iac.es/iac-star/ 
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errors in the derived luminosities. 

Translation to total stellar mass is more problematic 
due to a numerical issue in the lAC-STAR code that, 
for some sets of input parameters, gives a NaN result 
for some integrated quantities of the stellar population, 
such as extant stellar mass (A. Aparicio 2010, private 
communication). Using the stellar mass when available, 
and extrapolating those mass-to-light ratio trends with 
metallicity and age to those sets of parameters for which 
the stellar mass was unavailable, we find that the to- 
tal stellar mass per TRGB star is between 1.2 x 10'* 
and 2.1 x lO** Mq/Ntrgb, increasing as a function 
of both metallicity and age. We adopt a scaling of 
1.6 X 10^ Mq/Ntrgb, and again caution that different 
stellar populations could introduce systematic errors in 
the derived masses of order ~ 50%. 

5.2. Halo luminosity and mass 

The total number of TRGB stars we derive in our field 
is TVtrgb — 3.5 X 10*, resulting in a minimum halo lu- 
minosity of Lhaio > 4.6 X 10* Lq and minimum halo 
mass of Mhaio > 5.6 x 10* M©. We can estimate the 
total halo luminosity several ways. First, we note that 
we have covered approximately 30% of the area of the 
galaxy to a distance of 30 kpc, and that the halo den- 
sity drops sufficiently steeply that the majority of the 
stars should be contained within 30 kpc, resulting in a 
total of ihaio = 1.5 X 10^ Lq or Mhaio = 1.9 x 10^ Mq. 
Another estimate of the total luminosity can be deter- 
mined from the elliptical power law fit to the TRGB 
radial number density; however, as the slope of the pro- 
jected density is significantly more negative than —1, the 
integral diverges at small radius. We have therefore im- 
posed an inner cutoff at a radius along the minor axis 
of 5 kpc, and caution that this estimate does not ac- 
count for most of the halo stars that overlap with the 
galactic disk. The total halo luminosity using these as- 
sumptions is -Lhaio = 2.6 ± 0.3 X 10^ Lq, while the halo 
mass is Mhaio = 3.2±0.4xl0'' Mq. The main differ- 
ence between these two estimates is due to small-scale 
structure, rather than the counts beyond 30 kpc which 
only account for 15% of the luminosity using the power 
law fit. The quoted random errors, which are due to 
the uncertainty in the power law fit, certainly underes- 
timate the true uncertainty; for example, systematically 
shifting the inner points by the estimated systematic un- 
certainty of 0.2 dex, as in Section 14.11 changes the lu- 
minosity and mass by 50%. The confluence of these 
various estimates suggests that the total halo luminos- 
ity is -Lhaio ~ 2 ± 1 X 10^ Lq and the stellar mass is 
Mhaio ~ 2.5 ± 1.5 X 10^ Mq. 

We estimate the fraction of the total light and stellar 
mass of the galaxy contained in the halo from the total 
RC3 (Ide Vaucouleurs et aLlHool B and V magnitudes 
of B = 8.04, and V = 7.19, which gives an absolute 
magnitude of My = —20.6 and a F-band luminosity of 
« 1.5 X lO^*' Lq. This corresponds to a total stellar 
mass of « 4.4 x 10^° Mq, assuming M/Li 



riven by the B — V vs. M/Ly relation in lBell &: de Jona 
(|2001h . Therefore, a fraction ~ 0.13 of the luminosity 
and ~ 0.06 of the stellar mass of NGC 253 is in the form 
of its halo. This is consistent with theoretical models 
that predict that the median galaxy at this luminosity 
contains a stellar halo with 1 - 5% of the galaxy stellar 



mass (IPurcell et al.l[2007l) . 

The only other galaxies for which the total luminosity 
and mass of the stellar halo can be obtained are the Milky 
Way and M31, both of which are of similar luminosity to 
NGC 253. Measurements place the total luminosity and 
mass of the Milky Way's stellar halo at Ly ^ 10^ Lq 
an d M ^ 2 X 10^ Mq resp ectively (see the summary 
in lBullock fc Johnstonll200"5l) . an d the luminosity of the 
M31 stellar halo at L ~ 10^ Lq (llbata et al.ll2007| ). sim- 
ilar in magnitude to but somewhat smaller than that of 
NGC 253. 

6. CONCLUSIONS 

We have detected resolved RGB stars in the stellar 
halo of NGC 253 in both HST/ACS and ground-based 
Magellan/IMACS imaging data. The clean HST data are 
only able to probe a small area of the halo, but provide a 
benchmark by which we can evaluate the ground-based 
data. The large area covered by a single IMACS pointing 
provides a global picture of the halo, reaching out to 
distances of 30 kpc from the center of the galaxy, which 
is required given the substructure and large-scale features 
predicted in models of stellar halo formation and seen in 
the observations. 

While previous studies of resolved stellar halos have 
successfully detected stars in the halo and disc overed 
spectacular tidal features ()Mouhcine et al.l I2O1O0 . only 
for the Milky Way and M31 has there been significant 
quantitative analysis of the halo profile and structure. In 
this work, we have performed detailed artificial star tests 
in order to make quantitative measurements that can be 
compared to theoretical predictions and compared to ob- 
servations of other galaxies. 

The projected stellar density declines with elliptical 
radius as a power law of slope —2.8 ± 0.6. This is con- 
sistent with, but at the steep end, of the density profiles 
found for the M ilky Way and M31 halos; for example, 
iBell et all ()2008l) find that the three-dimensional density 
falls off as a power of index of between —2 and —4, cor- 
responding to projected radial power law indices of — 1 
to —3, and the projected radial profile of the halo of M31 
has a power law index of « — 2 (Jbata et al..2007 'l. An ex- 
ponential vertical distribution is also found to provide an 
acceptable fit, with an extraordinarily large scale height 
of hz = 5.3 ± 0.8 kpc. However, interpretation of the 
observed structure as a thick disk rather than a power 
law halo would require an unrealistically large vertical 
velocity dispersion, and is inconsistent with the observed 
halo density in the HST/GHOSTS fields. 

The halo is found to be flattened in the same direction 
as the stellar disk, with a projected axis ratio of « 0.35 ± 
0.1. This is consistent with what has been found for other 
galaxies: the Milky Way halo is also fla ttened in the 
same sense as its disk at s mall radii (e.g. iCarollo et al.l 
I2007[ ). lde Jong et all (|2009D find a typical halo axis ratio 
of c/a w 0.4 based on HST/GHOSTS data of nearly 
edge-on galaxies, and although not quantitative, the star 
counts in the halo of NGC 891 also appear s ignificantly 
flattened hkc the disk (jMouhcine et al.ll2010[ ). The dark 
matter halo that dominates the potential is also predicted 
in galaxy formation simulations to be flattened in the 
same s ense as the galact ic disk within 10% of the virial 
radius (jBailin et al.ll2005l ). Observations of more galaxies 
are required to see if this is a universal feature of stellar 
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haloS; 

Da vidgd ()2010[ ). who detected AGB stars in the halo 
of NGC 253, found their distribution to be flattened sim- 
ilarly to what is found in this study for the older RGB 
stars. Unfortunately, our observations are not sensitive 
to these younger stars, and so further observations will 
be required in order to fully understand the relationship 
between these populations. 

Density maps of the halo structure show substructure 
on a variety of scales. The shelf-like feature seen in earlier 
deep photographic studies (Malin & Hadley 1993) is evi- 
dent. In addition, there is a significant amount of smaller 
kpc-scale structure. These features are not apparent in 
the small HST fields but require the full IMACS field of 
view to appreciate. 

We estimate the total luminosity of the stellar halo 
to be ^ 2 ± 1 X 10^ Lq, and the stellar mass to be 
~ 2.5 ± 1.5 X 10^ Mq. This corresponds to a fraction 
0.13 of the total luminosity and 0.06 of the total stellar 
mass of the galaxy. These values are con sistent with the 
theoretical models of lPurceh etall (|2007[) . The NGC 253 
halo luminosity is a factor of w 2 larger than the halos of 
the Milky Way and M31, despite the similar total lumi- 
nosity of these galaxies, indicating significant galaxy-to- 
galaxy scatter in the relative halo fraction, as predicted 
by models. 

This paper provides a roadmap for analyzing ground- 
based imaging data of the stellar halos of galaxies beyond 
the Local Group in a quantitative way. These techniques 
will be applied in future work to a series of other galax- 
ies for which we have obtained similar data. The result- 
ing collection of quantitative global measurements of the 
scale and structure of stellar halos of a broad sample of 
galaxies will allow for statistical comparison to models of 
stellar halo formation, and will provide vital constraints 
on the galaxy assembly processes that are traced by the 
halos. 

A catalog of the detected sources is being made avail- 
able as a table on the Centre de Donnees astronomiques 
de Strasbourg (CDS). 
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